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ABSTRACT 

A Monte Carlo method is described for evaluation of the spectrum, 
directivity and polarization of x-rays diffusely reflected from 
stellar photospheres, The accuracy of the technique is evaluated 
through comparison with analytic re^^ults. Using the characteristics 
of the incident x-rays of the model for solar x-ray flares described 
in the accompanying paper, the spectrum, directivity and polariza- 
tion of the reflected and the total x-ray fluxes are evaluated. Tlie 
results are compared with observations. 



Also Department of Applied Physics 



I . INTRODUCTION 


In all models for Impulsive solar x-ray bursts the observed 
x-ray flux is partly duo to x-rays emitted directly toward the Earth 
and partly due to the reflection of x-rays emitted toward the photo- 
sphere. In models where the x-rays are produced by a beam of electrons 
directed toward the photosphere (see Petrosian, 1973^ paper I, and the 
accompanying paper III), the bulk of the x-rays (especially at high 
energies) are emitted toward the photosphere and the contribution of 
the reflected x-rays to the observed flux becomes significant. 

The general problem of the x-ray albedo of stellar atmospheres 
and its application x-ray binaries (treated in an approximate 
analytic manner by Basko, et will be discussed elsewhere. 

In this paper we present the expected characteristics ( spectrum, 
angular distribution and polarization) of the reflected and total 
(direct plus reflected) x-ray flux from the model described in the 
accompanying paper. The Monte Carlo procedure used in the calculation 
of the reflected x-rays is described in section II. Our results and 
their comparison with previous results (Tomblin, 1972; Santangelo et ai, 
1973; Henoux, 1975; Felsteiner and Opher, I976) are presented in 
section III. In section IV we summarize the results and comment on their 
comparison with the observations of the variation of the properties of 
impulsive x-ray properties with solar longitude. 

II. X-ray Albedo of the Solar Photosphere 

Ttie problem of diffuse reflection of x-rays from the photosphere 
is treated by using a Monte Carlo technique where the changes in the 
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properties of photons (with various initial energies, directions and polari- 
zations) are followed as the photons undergo Compton scattering in the 
photosphere before they escape or are absorbed by photoionization. 

Each photon is characterized by its momentum vector k = k(sln0 coS(p, 
sin 0 sin(p, cos 0 ), its polarization vector 8 (or the polarization angle 
between e and the vector Ic x z x It) and its depth z in the atmosphere. 
Instead of z it is more convenient to use the electron column density 
above z (n = electron density) 

N(z) = n(z' )dz' (1) 

25-2 

At the column density N 10 ''^cm where the photons undergo 
most of their scattering the density scale height is much smaller than 
the solar radius so the semi-infinite plane parallel atmosphere 
approximation will be used, and therefore the x and y coordinates 
of the photon need not enter into the problem (cf., however, the dis- 
cussion at the end of this section) . To further simplify the problem 
it is assumed that the radiation pattern of the x-ray source has 
azimuthal symmetry. Tills amoimts to assuming that the magnetic field 
lines (about which the electrons producing the x-rays spiral, cf. paper III) 
are radial, this removes any dependence on the azimuthal angle of the 
direction of observation. 

For the energy range (I 5 to 5 OO kev) of interest coherent scattering 
can be neglected ( cf. Viegele, Henry and Tracy, I966) and the scattering 
electrons (with binding energy w thermal energy in the range 10 to 100 ev) 
can be assumed to be at rest. Thus the Klein-Nishina formula can be 
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used to describe the scattering process 



(S 



2 2 2 
cos g cos cp + sin {p 


where C; k and , k' are the photon polarization vector and 
energy before and after the scattering through relative scattering angles 
B and tp. For energies below 'JO kov the classical Thomson cross section 
can be used for da without loss of accuracy. 

Tlie absorption process is described by a photoionization cross 
section which obeys a k power law between K-edges and has 

discontinuous jumps at the K-edges of the various elements. It is 
convenient to normalize all cross sections to the total classical 
Thomson cross section and measure column densities in units of the 
reciprocal of this number. The photo ionization cross section may then 
be written 




-a / k \-3 

i ^12.2 kev j 




k < k 


1+1 




(3) 


The values for are taken from Fireman (197^), with 01^ constant 

between the k-edges at k^ and At 12.2 keV Oi^ s 1.0 and the 

total effective photoionization cross section is equal to the total 
Compton scattering cross section. 

The Monte Carlo Procedure . Consider a photon after the i^ 
scattering at a column density with momentum and polarization 

vectors and characterized by energy k^ and angles 0^, 

♦ 

and relative to the fixed axes described above (for incident 

photons i = 0, Nq = O) . 

We first calculate the column density, 

^i+1 ”i “°®®i ' 

where the next event (absorption or scattering) occurs. The column 
density 6N traversed between events is obtained from the solution of 



where is a random number between zero and one. Two column densities 

^^kn ^%h generated using the total Compton scattering and 

photoionization cross section and respectively. The column 

density traversed is chosen as 6N = min(6N. ,6N and the process 

Kii ph 

associated with the smaller column density occurs unless N turns 

i+1 

out to be negative, in which case the photon escapes before the next 
event. If absorption occurs the position and energy deposited by the 
photon are recorded. 


k 


In the event of scattering the relative scattering angles Q and 
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are obtained from 




X ^0^'^kn^ ( 9 ^ ^2it) 'Sine' 

Jo (e' /2rt) *sine' 


( 6 ) 


and 


n = 


a,J(9,cp) 




(7) 


where R„ and R are random numbers between zero and one, and 
y cp ' 


^kn -J 


'P’ 


dO 


d(p 


(8) 


These integrals are inverted using a modified regulae falsi numerical 
method (Conte and DeDoor^ 1972). 'The energy of the scattered photon is 


k. 


i+1 




'm c 
e 


The polarization angle (|r relative to the coordinate with z-axis 
parallel to (i.e., the angle between and >< X 

is determined from an examination of the scattering cross section which can 


be written as 


da, 3c 


kn 


dn 


32if \ k / ^Pu Pl' * 


(9) 


with 


^u " k. 


idr 


'i+1 


(10) 
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The p portion of the cross section corresponds to unpolarized photons 
and is a relativistic correction. Tlie classical term, is a 

standard dipole scattering contribution. The polarization vector of the 
scattered photon, > must therefore lie in the plane of the old 

polarization vector and the new direction of propagation 

Therefore, if the photon is part of the p^^ (linearly polarized) 
component, \jr must satisfy 

Thus, p should be evaluated at this value of ijr before generating 
0 and cp ; 

Pj^ = 4( cos^0 cos^tp + sin^tp) , (12) 


The value of ilt is set by generating yet another random number, u , 

Pu ^ 

between zero and one. If R, > — r — ) then is chosen randomly 

'1 ' Pu'*'Pl 

between zero and 2n; otherwise i) is given by (11) . 

The values of the angles 0^, tp^ (relative to the fixed 

coordinates) after each scattering are kept track of by using a rotation 
matrix tp ) which is the ordinary rotation matrix for Euler 

angles 0^, {p^ (see e.g,, Mathews and Walker, I 97 O) and satisfies 


r^^,0,cp) = V^(ep) ^~\ q ) \ 

R"^(^, 9 ,cp)k = z (13) 

^ (')^e,;p)c s X . 
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The geometry of the system is shown in figure 1, 

A 

When a photon scatters the momentum vector polarization 

A 

vector given by 


A ^ 

= R (i|f,0,{p) 




The Euler angles of the scattered photon relative to the fixed xyz 
coordinate system, obtained from 




The Generation of the Incident Photons . The Incident photons are 
generated Unifonnly in the chosen energy band, uniformly in cos 0 q within 
the cone subtended by the photosphere at the x-ray production site, and 
uniformly in cp^ in the range 0 to 2n, Each photon is assigned 
a statistical weight equal to the photon number flux in direction 0 q 
at energy in the incident beam 

Vq = j(kQ) D(kQ,0Q)/kQ (16) 

where is the initial photon energy, J is the integrated flux and 

D is the directivity (see paper III for the definition of J and D 
and their form for the model of solar x-ray flares under consideration) . 
This weight is carried with the photon and is used to increment the 
proper bins in the histograms when the photon escapes. 
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Tho initial polarization angle should be goubrafcod so as to 

givo; on tho average; tho degree of polarization P( 9^) of tho 
incident radiation. In general; the incident radiation will bo partially 
linearly polarized and may be decomposed into a combination of Unpolnrlzod 
and linearly polarized light. Tlie linearly polarized and tho unpolarized 

components have relative weight [ P(ICq;9qJ j and 1 - | P(k^;0Q) | 

respectively. Tlie degree of polarization p(kQ;0Q) for the model of 
impulsive solar x-ray bursts under consideration hero is given in 
paper III. As described there the symmetry of the bremsstrahlung process 

A 

and the choice of the z - k^ plane ns the = 0 reference plane 

ensures that for P(1 <q;0q) <0 the linearly polarized component lies in tho 

direction = 0; and for P(kQ,0Q) > 0 it lies in the direction ([fQ = 

A A 

(normal to tho z - plane). Tlie polarization angle of the incident 
photon is picked by generating a random number; between 0 and one. 

If Ri)q < I P(kQ;0Q) I the photon is linearly polarized (with = 0 
for P < 0 or — for P > 0); otherwise is chosen randomly between 

0 and P.n, 

As mentioned at the beginning of this section the scattering has 
been treated in the plane parallel approximation. However; if tho source 
height above the scattering level is comparable to the radius tho proper 
transformation between coordinates referred to the surface nonaal at the 
flare site and those referred to the normal at tho point of entry into 
the photosphere should be made; and the inverse transformation should be 
made for escaping photons. 

This process is carried out by rotating the matrix; R('^q;^q; 9 q); 
of the photon in flare site Coordinates through the angle § between 
the normal at the flare site and the normal at the point 
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I 


1 


II 

of entry into the photO)?nhoro giving tho matrix n(^Q;0Q;O) to bo used 
for tho plane parallel approximation as 

jl 

n"\o,0,^Q) (17) 

A photon escaping after the nth scattering, ' must bo 

referred back to flare coordinates 

n<hy,'^) - R(V®n'«Pn^ J^(0;5,9o) 

Tho properties of each escaping photon are recorded in the proper 
histogram bins using the final direction; 'S, polarization; and 

energy; k ( the symmetry of the incident flux assures that there is no 
^ dependence) . 

III. RESULTS 

A. Comparison with Previous Work 

Before presenting the spectrum, directivity and polarization of the 
total (reflected plus direct) flux for the model of impulsive solar 
x-ray bursts described in paper III, we can test the Monte Carlo proce- 
dure and evaluate its accuracy through comparisons with earlier work. 

The problem of scattering from an atmosphere of electrons has been 
treated previously by several authors using both the classical and 
relativistic cross sections. Our results can therefore be compared to 
exact results from the theory of radiation transfer; approximate 
analytic solutions to the problem of a Compton scattering atmosphere 
and to previous Monte Carlo solutions. 
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Chanclrasokiiar (I96O) found the directional and polarization proper- 
ties of the z’adiation diffusely I'oflectod by a plane pai’allol atmosphere 
of perfect Rayleigh (dipole) scatterers, Abhyankar and Fymat (I97I) 
have extended this calculation to include scatterers with single 
scattering albedo = 0.99* 

Wo have found that with a few thousand pliotons the Monte Carlo 

procedure gives results in good agreement with the analytic solutions. 

An example is shown in figure 2 whore wo plot the directivity and 

degree of polarization of the radiation I’oflocted from an azimuthally 

symmeti'lc (j^ chosen randomly between 0 and 2rt) incident beam with 

0Q £= 60° and = 90°. The directivity is defined as the I’atio of 

the intensity in a given direction to the average intensity and the 

deg.‘'ee of pol^ ^i;.;iion is the difference between the intensity with 

polarization perpendicular to and in the meridian plane as a fraction 

of the total intensity; the mci'idiau plane colltalnD the dlrertion of 

observation and the noi’mal to the scattering surface. Here we have 

ignored the Compton shifty i.e.^ we have set a s Om ) and set 

101 X 

Oph = cr^/99 to give a single scattering albedo a?Q = O.99. The direc- 
tivity results are good to within the one sigma error bars inferred 
from the Monte Cax’lo results, and the degree of polar! ^jation results 
are good to much better than one sigma. 

Basko et al , by making various r.implifying assumptions, 

found an analytic approximation for the radiation diffusely reflected 
by a plane pai’allel Compton scattering atmosphere. Felsteiner and 
Opher (1976) have treated the same pi'oblem with a Monte Carlo approach. 
The spectrum of tJio radiation diffusely I'eflected in tlie backward 
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direction from a beam Of Tinpolarizod, normally Incident 30 photons 
is shown in fiffui’o 3. Our results agree very well with Felsteiner and 
Opher. Basko et al’s approximation III (shown in the figure) locates 
the single scattering peak and gives about the right height and width 
(Basko et al used a spread of incident energies between 29 and 30 keV) . 

We find an albedo consistent with Felsteiner and Opher and significantly 
less than Basko et al which is due to their lower iron abundance. 

Tlao problem of the radiation diffusely reflected by a Compton 
scattering atmosphere from a model spectrum (as opposed to a monochro- 
matic line) has been treated by Tomblin (1972); Santangelo et al (1973); 
Henoux ( 1975 ) and Felsteiner and Opher (I 976 ) . Ibmblin dealt with 
energies lower than we treat and Felsteiner and Opher used a thermal 
bremsstrahlung spectrum which is of interest in the case of x-ray 
binaries. Santangelo et al treated only the case of an isotropic power 
law source and did not report results for polarization. In figure 4 
we compare our results for the dix’ectivity with those of Santangelo et al 
and Henoux for roughly equivalent isotropic sources. At the lower 
energy the three curves lie within about 2 a of each other while at 
the higher energy our results agree very well with those of Henoux and 
also with Santangelo et al, taking account of the slightly different 
energy bands. 

B. Results for Solar X-ray Bursts 

Using the Monte Carlo procedure described in section II we have 
calculated the spectrum; directivity and polarization of the reflected 
radiation from the incident beam of the r;odel described in the accom- 
panying paper. Adding the reflected flux to the direct flux we obtain 
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tho spectrum, directivity and polarization of the total flux expected 
from this model. These predictions for the observable flare properties 
are presented below. 

i) Spectrum 

Tho spectrum of the total flux expected from a power law electron 

dJg .*5—1 

beam (with a differential number flux « e“°” , c.f. paper I, eq. l*f) 

is presented in figure 5. In order to show the detailed properties of 

the spectrum we have plotted n(k) versus photon energy where n(k) = J(k)/k 

is the photon number flux (per sec per unit energy interval) . As is 

evident from these curves the steepening of the spectra beyond 100 keV 

which was present in the direct flux (c.f. papers I and III) becomes 

more pronounced when tho reflected flux is taken into account. Such 

breaks have been observed (Frost I969 and Frost and Dennis 197^) during 

the impulsive (but not the slow) phase of tho flare for flares with 6^3. 

Our model shows that a larger change in the spectral index occurs at 
100 keV for harder flares than for softer ones, although experimentally 
the determination of a break for softer flares would be more difficult. 

In Table I we present results from least square fits of the spectra 
to two power laws, one for photon energies between 20 and 100 keV and 
another for photon energies greater than 100 keV. The residual variance 
is much less for the two power law fit than for a single power law fit 
across the entire energy range. Using the F-test it is possible to test 
the hypothesis that the data is actually a sample from a single power 
law with statistical fluctuations and that the reduction in variance 
with the two power law fit is due simply to the increase in the number 
of free parameters in the fit. The results show that for none of the 
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curves can the achieved reduction be expected more than of the time 
on a statistical basis. Taken as a whole this strongly indicates, but 
by no means proves, that our Jnodel spectra are best fit by a two power 
law spectrum. 

Sometimes such deviations of flare spectra from a single power law 
are interpreted as indications of a thermal bremsstrahlung of the flare 
x-rays (c.f., e.g. Elcan I 975 ) . We have therefore fitted our model 
spectra to a thermal spectrum and determined the temperatures shown in 
Table I , Tlie thermal spectrum provides a poorer but nevertheless 
acceptable fit. 

We would also like to point out the systematic flattening (or 
hardening) of the spectrum (in 20 to 100 keV range) of harder flares 
(6 S 4) in going from limb flares (cos 0 = 0 . 0 ) to flares at the 
center of the solar disk (coS0 = -1.0). Because the majority of 
flares have 6 between 3 and 4 , this model indicates that on the 
average limb flares would be steeper (or softer) than those at the 
center of the disk although the variation is not as strong as observed 
(c.f, however, section IV) . 


TABLE I 


Results of Fits to the Spectra 




|i S COS0 


0 

0 s: U 5 S -0.2 

- 0.4 s se -0.6 

-0,8 s |i s -1.0 

* 

X 

3 

.55 ± .03 

.40 ± .06 

.20 ± .08 

20 < k < 100 keV 

4 

.57 ± .02 

.61 ± .08 

.46 ± .09 


5 

.6l ± ,01 

.69 ± ,06 

.61 ± .11 

* 

X 

3 

1.64 ± .04 

2.68 ± .26 

3.5^ ± . 9 ^ 

k > 100 keV 

4 

1.52 ± .03 

2.25 ± .07 

2.86 db .54 


5 

1.52 ± .01 

1.98 ± .07 

2.35 ± .16 


3 

20.7 ± 3 

20.3 ± 3 

19.1 ± 3 

kT keV 

4 

16.4 ±3 

15.9 d; 3 

15.7 ± 3 


5 

13.6 ± 4 

13.1 ± 4 

12.6 ± 4 


■X- 

X 


d loe J(k) , , . 

d log ir" +5-1 


j(k) = photon flux in 


keV 


cm -sec-keV 


Error bars are for a 99 ^ confidence interval. 























ii) Directivity 

The directivity D of the total outward (-1 < cos0 < O) flux 
normalized at cos0 = -1 (center-of -disk flare) is plotted for various 

fsf 

energies and index 6 in figure 6. (Note that D{k;0) =f 
is defined differently than the directivity D in paper III.) Compai’i- 
son of these curves with the directivity of the direct flux of paper III 
shew that between 20 and lOO keV the radiation is usually fairly 
isotropic. For example, at 20 keV the ratio of center to limb flux 
is 1:6 for the direct f 1' jc while for the total flux, in the case of 
zero pitch angle and an electron energy spectrum with index 6 ” 3; ^ 
and 5; these ratios are 1:1, 1:1.5 ^nd 1:1.8, respectively. Tliis 
structure is even moi’e pronounced at 100 keV, and for 30 < k < 100 keV 
and 6=3 the limb flares (cosg = O) become weaker than the centor-of- 
disk flares. 

Above 100 keV the scattered flux makes a smaller relative contri- 
bution due to larger relative energy shifts of the scattered radiation 
up the steep power law spectrum. The scattered flux also makes a 
reduced contribution below 20 keV due to stronger absorbtion through 
photoionization. Tlie result is that outside 20 to 100 keV the direc- 
tivity is similar to that of the direct flux. Tlie directivity does 
not depend too strongly on the pitch angle but does tend to decrease 
with Increasing pitch angle. 

iii) Polarization 

The degree of polarization of the total flux is shown in figure 7- 
In general, as comparison of these curves with the curves in paper III 
for the degree of polarization of the direct flux shows, the scattered 
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flux reduces the degree of polarization from that of the direct flux, 
but the overall shape Is preserved. The degree of polarization 
decreases both with increasing pitch angle and increasing energy 
but Increases with Increasing 6, the power law index. The scattered 
radiation affects the degree of polarization only slightly because for 
cosi? w -1, where it dominates, both the direct and scattered flux have 
zero polarization, while at cos0 0, where the direct flux is 
strongly polarized, there is little scattered flux, 

iv) Albedo 

The albedo of the solar atmosphere for a monochromatic incident 
flux is shown in figure 3* In general, for steep spectra of incident 
flux, such as those occurring in solar x-ray bursts, at a given energy 
most of the contribution to reflected flux come from incident 
photons of similar and slightly higher energies with little contribution 
from much higher energies. Thus, in general, monochromatic albedos 
will be good approximation for steep power law incident flares. Because 
of strong directivity of the incident flux of a solar flare, albedo 
definitions become complicated. On table II we give D _(k,ii) the 
reflect photon number flux as a fraction of the total (integrated 
overall energies and directions ,yj(k)dk) generated flux. For comparison 
we also show this fraction for the direct flux which is the same as 

rw 

D(k,cos0) [note that D + oc D] . As evident for center-of-dlsk 

flares the reflected flux is 5 to 6 times stronger than the direct 
flux while at the limb the reflected flux contributes only about 
35^ to the total observed flux. 
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The observed flux for a given incident electron beam is given by 

m 

= [D(k^0) + J(^) • (19) 

Thus with the help of this equation and the relation between J(k) 
and the incident electron beam one can calculate the total energy of 
the accelerated electron for observed fluxes at any value of 0, 

We also note that most of the photons not reflected are absorbed 
within three scattering lengths. 
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TABLE II 


Comparison of Direct and fieflected Fluxes at 20 - 25 keV 



6 

" 1 

0 < p < -0.2 

-0.8 < li < -1.0 


3 

.o 4 o 

0.122 

nj 

.052 

0.153 

5 

.054 

0.129 

D(k,p) 

3 

.078 

0.024 

B 

.091 

0.022 

5 

.094 

0.021 


p. = COS0 


/ ' -2 ~1 

J(k)dk = total incident flux (keV cm sec ) 
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IV. SUMMARY AND CONCLUSIONS 

We have developed a Monte Carlo procedure for evaluation of the 
spectrum, directivity and polarization of x-rays reflected from 
an attnosphoro containing electrons and ions. Wo have tested 
the accuracy of this procedure by comparing its results with analytic 
solutions of similar problems. Our results also agree with previous 
numerical and Monte Carlo solutions. Tliis program has been used to 
determine the properties of the reflected x-rays from the solar atmosphere 
during a solar flare. Using the polarization, spectrum and directivity 
of the x-rays expected from the model described in paper I and the 
accompanying paper (paper III), we have calculated the expected 
properties of the total (reflected plus direct) radiation for various 
parameters of the model. In comparing the total flux \vith the direct 
flux (the flux in absence of reflection) we find the following: 

1) The total flux is affected mainly in the 20 to 100 keV range 
of energies. Above 100 keV the contribution of the reflected radiation 
is small because of increasing Compton shift with energy, and there is 
little reflected radiation below 20 keV because of strong photo- 
absorption at these energies. 

2) Because of the above the steepening of the spectrum above 
100 keV is more pronounced than that reported on paper I (and is in 
better agreement with observations by Frost (I 969 ) for the direct flux). 

We also note that lor flares with smaller electron index 6 (i.e. harder flares) 
the spectra tend to be steeper near the solar limb, cos0 = 0, in comparison with 
those at the center of the disk, cos0 = -1. This phenomenon disappears 
as 6 increases. In some cases the curvature produced by this break 
could result in a spectrum (from a power law electron spectrum) resembling 
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a thermal bremsstrahlung spectrum Indicating that observation of such i * 

spectra does not necessarily indicate a thermal origin of the radiation, 

3) nie directivity of the total flux in the 20 to 100 koV range ^ 

is reduced considerably in comparison with the directivity of the 

direct flux. For electron spectra with index 6 « 3 the flux is isotropic 

at 20 - 30 keV; and has the opposite sense of directivity (stronger in the i 

center than at the limb) compared to the direct flux in the 30 - 1^0 koV 

range of energies. For higher values of 6 directivity is similar to 

that of the direct flux but the ratios of limb to center fluxes are 

smaller. 

4) In general^ the total flux is still highly polarized but the 
degree of polarization is smaller than that of the direct flux. 

In paper II (Petrosian 19T^) the properties of the direct flux 
were compared with the statistical behavior of a sample of flares ' 

described by Datlowe, Elcan and Hudson (197^) ■ Two important aspects 
of this comparison were the absence of variation of frequency of 
occurrence of flares across the solar disk and steepening of the spectra 
toward the limb. It was shown that although for the direct flux the 
spectra are identical across the disk the strong directivity of the 
direct flux causes favorable conditions for observation of steeper 
spectra toward the limb. As mentioned above the directivity of the 
total flux is smaller or non-existent at the 20 - 30 keV range. Howevei*, 
as mentioned in items (2) and (3) above, lor values of 6 where directivity 
is small the specti'a of the total flux tend to be steeper (or softer) on 
the limb, while for higher fl where this steepening is absent the flare 
tend to be stronger toward the limb than at the center. 
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The situation hoi*o is more complicated than that of the direct 
flux, and the simplified analysis of paper II becomes complicated and 
not worth being repeated because of its simplicity. However, it can 
be shown that the above two effects tend to result in gradual steepening 
of the spectra toward the limb. In fact, an analysis similar to that 
of paper II shows that for the total flux one can define an equivalent 
parameter b so that the b/a 6 ^ term in paper II is about unity 
implying a difference in spectral index between flares at -1 < cos 0 < -|- 
and at < cos 0 <0 of about +O .5 Instead of the observed 0,7. Note 
also that because of the smaller value of b/afi' the expected fre- 
quency of occurrence of flares would be more uniform across the disk 
(once the visibility factor is taken into account) than for the 

direct flux. 

The negative polarization is in agreement with observations of 
Tindo et al (1971); but because of the tentative nature of these 
observations not much more could be inferred from their comparison 
with the model. 
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FIGURE CAPTIONS 


Figure 1 Tlio Eulor angles tpf and ij; which transform the x-axis 

into the photon polarization vector S; and the z-axis into 

A 

the momentum vector Ic. Adopted from Mathews and Walker (1970) 
Figure 2 Comparison of the Monte Carlo results, represented by circles, 
with the analytic calculations of Abhyanknr and Fymat (1970) 
(linos), for directivity (open circles and dashed lino) 
and polarization (filled circles and solid lino) of scattered 
radiation from a plane parallel atmosphere consisting of 
Rayleigh scatterors with single scattering albedo Wq = 0*99* 
The incident beam is normal to the atmosphere (cosG^ 0, 

uniformly distributed between 0 and 2jt) and hac polarize^ 
tion vector normal to the meridian plane, = 0, 
Appi’oximately 4,000 emerging photons. 

Figure 3 Albedo (normally reflected photon number per keV per sr 

per incident photon) versus photon energy of an atmosphere 
containing electrons and ions with "cosmic" abundance 
ratios for a normal incident monoenergetic beam of 30 keV. 

Tlie solid line is fitted to the filled circles from this 
paper. Tlie open circle from Felstoiner and OpUer (197^) 

Monte Carlo results. The dashed line from semi-analytic 
approximation III of Basko et ul (1974) with incident photon 
energies between 29 and 30 keV and for lower iron abundance. 
Figure 4 Comparison of the present Monte Carlo results (solid lines 
and points with error bars) with those of Henoux(l975) 

(dashed lines) and Santangelo et al (1973) (dotted linos). 
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All curves are for nn unpolnrized incident plioton flux of 


index 4 = 3 . Lower curves for reflected photon onorgios 
50 < k < 70 koV; middle curves for 80 < k < 100 koV; upper 
curves for I5 < k < 20 keV except for the clotted line which 
is for 15 < k < 30 koV photons. Santangelo et al results 
are for a lower iron abundance and for 6 = 2.5. 

Figure 5 Spectrum of total (direct plus reflected) photon flux for 

various index 6 of the incident flux. n(k) = photons por 
keV per sec; k = photon energy; the curves marked -1 give the 
flux for -1 < COS0 < - 0.8 (center-of-disk flares) and the 
curved marked zero give the flux for -0,2 < cos 0 < 0 (limb 
flares). Pitch angle T] = 0 , 

Figure 6 Directivity, normalized at cos 0 = -1 (center-of-disk flares), 

of the total flux for various index 6 and energies k in koV. Pitch 
angle 7 ] = 0 for all curves. For 7 ] 0 directivity is 

slightly smaller than for 7] = 0. 

Figure 7 Degree of polarization of the total flux versus cos0 for 

various values of the index 6, pitch angle 7] and energy k. 

In (a) all curves are for 6=4 and are labelled by the 
value of T] in degrees. In (b) curves are labelled by the 
value of index 6. The dotted line is for an unpolarized 
Incident (or direct) flex and the dashed line gives the 
polarization of the direct flux at 21 , keV. Points are 
from the observations of Tindo et al (I972) , 
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